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Abstract. In this paper we study the quasistatic crack growth for a cohesive 
zone model. We assume that the crack path is prescribed and we study the time 
evolution of the crack in the framework of the variational theory of rate-independent 
processes. 



I. INTRODUCTION 

In this paper we present a variational model for quasistatic crack growth in the 
presence of a cohesive force exerted between the lips of the crack. 

The evolution of the crack is governed by an energy which is the sum of three terms: 
the bulk energy of the uncracked part, the energy dissipated in the fracture process, 
and the work of the external loads. The main mathematical difficulty is given by the 
fact that the fracture energy depends on the opening of the crack. For this reason 
we cannot apply directly the tools developed so far in the applications to fracture 
mechanics of the theory of free discontinuity problems (see [10], [Hj, [Hj, P] [Hj, [Sj, 
H)- 

To simplify the mathematical difficulties, we assume that the crack path is pre- 
scribed, and we focus only on the time evolution. This allows us to consider very 
general bulk and crack energies, which may include constraints on the crack open- 
ing, related to the infinitesimal noninterpenetration of matter. The evolution of the 
crack is defined (see Definition 13.41 below) in the framework of Mielke's approach to 
a variational theory of rate-independent processes (see [T2], [TT] ). 

We prove an existence result for the quasistatic evolution, by approximating the 
continuous-time problem by discrete-time problems, for which the evolution is defined 
by solving incremental minimum problems. The irreversibility of the crack process 
leads to introduce an auxiliary time-dependent function t i— > j(t) (see Section El 
below), defined on the prescribed crack path, which takes into account the local 
history of the crack up to time t. The main mathematical difficulty in the proof 
is the compactness of the approximating functions t t— > 7fc(t). This is solved by 
introducing a new notion of convergence of functions related to the problem, with 
good compactness and semicontinuity properties. 

2. SETTING 

The reference configuration is a bounded open set Q of R n with Lipschitz boundary 
dfl, which can be written as the union of two disjoint Borel sets 8qQ and d\Q, with 
Jrff n ~ 1 (doQ) > and d±Q relatively open. Here and henceforth ffl n ~ x denotes the 
in — l)-dimensional Hausdorff measure. On doQ, the Dirichlet part of the boundary, 
we will assign the boundary deformation, while on d\Q, the Neumann part of the 
boundary, we will prescribe surface forces. 

We assume that the cracks are contained in a compact C 1 -orientable (n — 1)- 
dimensional manifold M C £1 with boundary dM, such that Q \ M is connected. 

l 
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Therefore it is reasonable to take the deformation u as a function in the space W 1,P (Q\ 
Af ; W 71 ), so that the essential discontinuity points of u are contained in Af. Although 
the natural choice is m = n, there are no mathematical difficulties in considering an 
arbitrary m > 1. The case m = 1 is used in the study of antiplane shears. The 
number p > 1 depends on the bounds on the energy density considered below. 

We take into account prescribed time-dependent boundary deformations t i— > ijj(t), 
with ip(t) G W 1,P (Q; IR m ), in the sense that for each time t G [0,T] we consider only 
deformations u G W l ' p {Tt \ Af; R m ) such that 

u = ip(t) on dofl, 

where the previous equality has to be considered in the sense of traces. We assume 
also that, as a function of time, t \— > ip{t) is absolutely continuous from [0, T] into 
W l ' p (n;R m ). 

Thus the time derivative 1 1-> ip(t) belongs to the space L\[0, T\; W x *(to', M m )) and 
its spatial gradient t i-> V^(t) belongs to the space £ x ([0, T]; L p (fi; M mxri )). 

We assume that the uncracked part of the body is hyperelastic and that its bulk 
energy relative to the deformation u G W 1,P (VL \ Af; W 71 ) is of the form 



W(x, Vu) dx, 

where £) is a given Caratheodory function W: (£1 \ M) x M mxn — > R such that 

(Wi) £ i— > VK(x, £) is quasiconvex and C 1 for every \ M; 

(W2) there are two positive constants a , ai and two nonnegative functions 60, &i G 
L l {n s Af) such that 

a ier - &o(x) < < ai ier + m*), (2.1) 

for every 6(H\ Af) x M mxn . 

Since £ i— » £) is rank-one convex on M mxra for every Af, we can deduce 

from (|2.1j) an estimate for the partial gradient of W with respect to £, c^W 7 : (fi\Af) x 
jypixn _^ jyjmxn^ jy[ ore precisely, there are a positive constant a 2 and a nonnegative 
function b 2 G \ Af) such that 

<02|er _1 + & 2 (x), (2.2) 

for every (x,£) G (fi \ Af ) x M mxn . 

To shorten the notation we introduce the function W : L P (Q \ M;M mxn ) — ► R 
defined by 



:= / W(x,*)dx, 

for every * G L p (Vl\M; M mxn ). By flTIJ) and ((221) the functional W is of class C 1 on 
L p (Q \ Af; M mxn ) and its differential dW: L P (Q \ Af; M mxn ) L 9 (ft \ Af; M mxn ), 
+ = 1, is given by 



for every $, \I/ G L p (fi \ Af; M mxn ), where (•, •) denotes the duality pairing between 
the spaces L 9 (ft s Af; M mxn ) and L P (Q \ Af; M mxn ), and d^W(x, *) : $ denotes the 
scalar product between the two matrices d^W(x, \I/) and $. 
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By the assumptions on W, the functions W and dW satisfy the following properties: 
there are two positive constants «o, a± and two nonnegative constants /3q, fti such that 

« 11*11; - < nw < «! ii*n; + a, (2.3) 

for every \l/ G L p (f2 \ M; M mxn ), and there is a positive constant «2 such that 

{m(*),$) <a2(l+\\*\\F X )\\*\\p> ( 2 - 4 ) 

for every $ G L p (fi \ M; M mx "). 

For a fixed time t G [0, T], we assume that the external time-dependent loads Jz?(i) 
belong to \ M;M m ))', the dual space of W 1 *^ \ M;R m ). The duality 

product (j£f(t), it) is interpreted as the work done by the loads on the deformation u. 

Let us fix an orientation of M and let it® be the trace of it on the positive side of 
M, and it® be the trace of u on the negative side of M. The most general form of 
the work done by the external loads is given by 

(jSf (t),u)= f f{t)udx+ I H(t):Vudx + 

+ ! g(t)udJ^ n - 1 + ( {g®(t)u® +g e (t)u e )dM> n - 1 , 

where f(t) G L 9 (fi \ M;R m ), H(t) G L 9 (fi \ M;M mxn ), g(t) G L 9 (difi;R m ), #®(t) 
and g e (t) G L 9 (M;M m ), with p" 1 + g _1 = 1. Actually the representation theorem 
for (W l,p (Q \ M;M m ))' shows that it is enough to use just the terms of the first 
line of (|2.5|) . The terms in the second line have been added in order to write in an 
explicit way the contribution of the surface forces acting on the Neumann part of the 
boundary and on one or both sides of M. 

With these assumptions we do not exclude the possibility that H(t) could be dis- 
continuous on M. Moreover, observe that if / (t) , H (t) , g(t) , g®{t) and g e (t) are suf- 
ficiently regular, then 

f(t) - div H(t) 
plays the role of the volume forces on Q \ M, 

g(t)+H(t)u 
plays the role of the surface forces on diQ, and 

g @ (t) - H®{t)v and g e (t) + H e (t)v 

play the role of the surface forces acting on the positive (respectively negative) side of 
M, where v is the outer unit normal to d(Q \ M). We observe that, by our positions, 
v turns out to be the inner normal on the positive side of M; this is why in the last 
formula we take the minus sign in front of H®(t)v. 

We assume that, as a function of time, t i— » Jz?(t) is absolutely continuous from 
[0,T] into (W^itt \ M;R m ))'. Thus the time derivative t i-» Sf(t) belongs to the 
space L\[0, T); (W llP (fi\M; R m ))'). If Jgf(t) is represented by (f2~K|) . then the absolute 
continuity of t (— > j£f (i) follows from the absolute continuity of the functions 1 1— > f(t), 
t ^ H(t), t ^ #(t), t h+ #©(£), and t h+ # e (t). 

If the deformation it has a nonzero jump [it] = it® — it® on M, then the body has 
a crack on (part of) M. More precisely the crack is given by the set 

{x G M : [u](x) + 0}. 
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Let us consider now the work done to produce a crack. If we neglect for a moment 
the problem of irreversibility, we may assume that this work can be written in the 
form 



' M 

where tp: M x M m — > [0, +00] satisfies the following properties 

(ipi) (p is a Borel function; 

(<p 2 ) (p(x, 0) = for Jg^-a.e. x G M; 

((ps) the function y 1— > <p(x, y) is lower semicontinuous on R m for t ^ 7n ~ 1 -a.e. x G M. 
A simple example is given by the function 

v(x ,y) := + «» E,V » (2.6) 



where a > and 6 > are real constants. The constant a plays the role of an 
activation energy; if b > 0, there is also an energy term proportional to the amplitude 
of the crack opening. The classical Griffith's model corresponds to the case a > 
and b = 0. 

Let L°(M) be the set of extended real valued measurable functions on M and let 
L°(M)+ be the set of functions w G L°(M) such that w > JT n_1 -a.e. on M. 
We introduce the function 0: L p (M;lR m ) -> L°(M)+ defined by 

0(u>)(x) := <p(x,w(x)), 

for every w G L P (M; R m ) and for n ~ x -a.e. x G M. 

Given an arbitrary family {wi)i & i in L°(M) + the essential supremum 

w = ess sup Wi 
iei 

of the family is defined as the unique (up to ^equivalence) function in L°(M) + 
such that 

• w > w i J^ n_1 -a.e. on M for all i G /; 

• if z G L°(M)+ and 2 > JT^-a.e. on M, then z > w Jf^-a.e. on M. 

For the existence of such a function see, for instance, [T^ Proposition VI- 1-1]. 

Suppose now that the deformation u depends on time, i.e., we have a map t h— > u(t) 
from [0, T] into W l,p (Q \ M; M m ). If no crack is present until time and 

<KKs)D < <t>([u(t)}) ^ n-1 -a.e. on M 

for every s G [0, t], then the energy dissipated in the crack process in the time interval 
[0, t] is given, in our model, by 



M 

This happens for instance when s i— > <p([u(s)}) is monotonically increasing ^ n_1 -a.e. 
on M. 

In the general case, the irreversibility of the fracture process leads to introduce an 
auxiliary function t \— > /3(t) from [0, T] to L 1 (M), which takes into account the history 
of the system up to time t. We assume that for every < t\ < t 2 < T we have 

P(t 2 ) =/3(ti) Vesssup^([u(s)]) ^ n_1 -a.e. on M, (2.7) 

tl<S<t 2 
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so that 

P(t 2 ) - (3(h) = ess su P (0([m(s)]) - (3(h)) + JT^-a.e. on M, 

t!<S<t 2 

where for every o G 1, fl + := a V denotes the positive part of a. 
In particular 

• t i-> f3(t) is increasing, i.e., /3(h) < (3(t 2 ) n_1 -a.e. on M for < t x < t 2 < 7; 

• P(t) > i([u(t)}) JT^-a.e. on M for every t G [0,T]. 

In our model the energy dissipated in the time interval [ti,^] is given by 
\\P(h) - (3(ti)\\i,M := [ ((3(h)- (3(h)) d^ 1 . 

J M 

According to this assumption there is no dissipation in the intervals 1 2 ] where 
0([w(s)]) < (3(ti) J^ n_1 -a.e. on M for every s G [ti, t 2 ] , while the dissipation is given 
by 

/ (<j>([u(t 2 )}) - ^(h)})) d^' 1 

J M 

whenever (3(ti) < (fi([u(s)}) < 4>([u(t 2 )]) for every s G [ti,^]- 

It follows from ()2.7|1 that (3(t) is uniquely determined by (3(0) and by the history of 
the deformation s \— > u(s) in the interval [0,t\. Since it is difficult to deal with ()2.7j) 
directly, we prefer to define the notion of quasistatic evolution by considering a more 
general internal variable t i— > y(t) which is assumed to satisfy the following weaker 
conditions: 

• t i— > y(t) is increasing, i.e., 7(^1) < 7(^2) J4? n -&.e. on M for < t\ < t 2 < 7; 

• l(t) > i([u(t)}) JT^-a.e. on M for every t G [0,T]. 

We do not assume from the beginning that t 1— > y(t) satisfies ()2.7|) . This property 
will be a nontrivial consequence of the other conditions considered in the definition 
of quasistatic evolution (see Theorem 13 .7|) . 

Given functions if) G W 1,P (Q; M m ) and 7 G L°(M) + , it is convenient to introduce the 
set AD(ij;,j) of admissible deformations with boundary value ip on <9ofi and internal 
variable 7. It is defined by 

AD(ip, 7) := {m G W^' p (fi \ M; M m ) :</>([«])< 7 on M, and u = ^ on <9 fi}, 

where equalities and inequalities are considered J4? n ~ l -&.e., and the last equality refers 
to the traces of u and ip on <9o^- 

An admissible configuration with boundary value ip on <9of2 is a pair (w,7), with 
7 G L 1 (M)+ := L X (M) n L°(M) + and w G AD(>,7). 

3. Definition and properties of quasistatic evolutions 

For every t G [0,T], the total energy of an admissible configuration (^,7) at time t 
is defined as 

<f(f)(w, 7 ) := W(Vw) - (J?(t),u) + || 7 ||i,m, 
where || • denotes the L 1 -norm on M. 

We now introduce the following definition in the spirit of Griffith's original theory 
on the crack propagation. 

Definition 3.1. A pair (u, 7) G W l > p (tt \ M; R m ) x L l (M) + is globally stable at time 
t G [0,7] if ue AD(ip(t),y) and 

7) <*(*)M) (3.1) 
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for every 5 > 7 and for every v G AD(i[)(t), S). 

In other words, the total energy of (u, 7) at time £ cannot be reduced by increasing 
the internal variable 7 or by choosing a new admissible deformation with the same 
boundary condition. 

Remark 3.2. For every £ G [0,7] let (u(t),~/(t)) G W 1 *^ \ M;R m ) x L 1 (M)+ be 
globally stable at time £. By Definition Kill we can deduce an a priori estimate on 
«(£). Indeed, by comparing <g(t)(u(t), j(t)) with <g(t)(ip(t), 7(£)), which is bounded 
uniformly with respect to t, we get that W(Vw(t)) — (J? (t) , u(t)) is bounded uniformly 
in time. Next, by the assumption ()2.3|) on W and the boundedness of jSf(t) in 
(W 1,p (fi s M;R m ))', we obtain that the W^-norm of «(£), ||«(t)||i,p, is bounded 
uniformly with respect to £. Furthermore from this fact and by Definition 13 . II we get 
that the crack term ||7(£)||i,a/ is bounded uniformly in time, too. 

Remark 3.3. Condition (|3.1jl is equivalent to 

£(t)(u, 7) <£(t)(v, T v 0(H)), 

for every v G W l > p (tt \ M;R m ) such that v = ip(t) Jif n ~ l -a.e. on doSl. This is 
equivalent to 

W(Vti) - <J2f(t),t*) < W(V«) - (J2f + ||(0(H) -7) + ||i,m (3.2) 

for every u G \ M; M m ) such that v = ip(t) Jf^-a.e. on d tt. This implies 

that if (u, 7) is globally stable at time £ and 7 G L 1 (M) + satisfies 0([w]) < 7 < 7 
J^ n_1 -a.e. on M , then (w, 7) is globally stable at time £. 

Definition 3.4. An irreversible quasistatic evolution of minimum energy configura- 
tions is a function £ ^ (u(t),j(t)) from [0, T] into W 1 *^ \ M; M m ) x L^M) 4 " which 
satisfies the following conditions: 

(a) global stability: for every £ G [0, T] the pair («(£), 7(£)) is globally stable at 
time £; 

(b) irreversibility: 7(3) < 7(£) n_1 -a.e. on M for every < s < £ < T; 

(c) energy balance: the function £ 1— > <f (£)(«(£), 7(£)) is absolutely continuous on 
[0,T] and 

|(<r(£)( w (£),7(£))) = (0W(Vu(f)), V^(£)) - (J2f (*),#)) - (^(£),n(£)) , 

for a.e. £ G [0,T]. 

Remark 3.5. Condition (c) is equivalent to the following one: 

(c') energy balance in integral form: the function £ 1— > (9W(V , u(£)), VV ; (£)) — 
(J2f (£),«(£)) belongs to ^([O.T]) and 

^(£)K£),7(£))-^(0)K0), 7 (0)) = 

= jf (<0W(Vu(s)), V^(3)) - (JSf (s), ^(s)> - (jSf (s), «(*))) 

for every £ G [0, T]. 
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This can be written in the form 

W(Vu(t)) - W(Vu(0)) + || 7 (t) - 7 (0)||i,m = 

= J ({dW(Vu(s)), V^(s)) - (JSf rf S + (3.3) 

+ (JSf(t), u(t)> - (jSf(O), u(0)> - / <JSf(s), u(s)> tte, 

./o 

for every t G [0, T]. The first line is the increment in stored energy plus a term which 
will be interpreted as the energy dissipated by the crack process in the time interval 
[0,t], as we shall see in Remark 13.81 Using the divergence theorem we can show that 
the second line represents the work done in the same time interval by the forces which 
act on 8qQ to produce the imposed deformation. The third line represents the work 
done by the imposed forces in the interval [0, t]; this follows from an integration by 
parts when t i— > u(t) is regular enough, and can be obtained by approximation in the 
other cases. 

If t \— ► («(t),7(£)) satisfies condition (a), then (u(t),j(t)) is bounded in W 1,P (Q \ 
M; IR m ) x L 1 (M) + by Remark 13.21 Therefore in condition (c') it is enough to assume 
that t i — ^ (<7W(Vu(t)), W>(t)) - (jSf(*),u(*)) is measurable. 

In the following theorem we prove one inequality of the energy balance. 

Theorem 3.6. Let t ^ (u(*), 7(f)) be a function from [0, T] into W X *(Q \ M;R m ) x 
L 1 (M) + which satisfies the global stability condition (a) and the irreversibility con- 
dition (b) of Definition]^ Assume that t ^ (dW(Vu(t)),Vip(t)) - {Sf(t),u(t)) is 
measurable. Then 

^(t)(«(t),7(*))-^(0)(«(0),7(0))> 
> J (<0W(Vu(s)), W(s)) - {&{s)A{s)) - (J2f(s),«(a)>) ds 
for every t G [0, T}. 

Proof. We note that t ^ (d\V(Vu{t)),Vip{t)) - (5?(t),u(t)) belongs to ^([0,^]) by 
the estimates of Remark 13.21 The result can now be obtained arguing as in jl] (see 
the proof of Lemma 7.1 and the final part of the proof of Theorem 3.15). □ 

Now we prove that for a quasistatic evolution t > (u(t),^(t)), the internal variable 
1 1 — ► 'y(t) satisfies a condition analogous to (|2.7|) . 

Theorem 3.7. Let t 1— > (u(t),j(t)) be a quasistatic evolution. Then 

7 (t 2 ) = T (t x ) Vesssup0([w(s)]) Jf n_1 -o.e. on M, (3.4) 

/or every < £1 < £2 < 7\ 
Proof. It is enough to prove that 

7 (t) = 7 (0) V ess sup <f>([u(s)]) jr n_1 -a.e. on M, (3.5) 

0<s<i 

for every £ G [0, T]. Let 7 (£) be the right-hand side of ()3.5j) . Since £ 1— > 7 (£) is 
increasing and </>([«(£)]) < 7 (t) Jf n-1 -a.e. on M for every t G [0, T], it follows that 
lif) < 7 (£) J^^-a.e. on M for every £ G [0,T]. As 0([n(£)]) < 7(£) jT^-a.e. on 
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M, by Remark E31 the pair (u(t), j(t)) is globally stable at time t for every t G [0, T\. 
Since t i— > 7(£) is increasing, we can apply Theorem 13.61 and we obtain 



for every t G [0,T]. By the energy balance (c) it follows that &(t)(u(t),'y(t)) > 



W(V«(t)) - (J2f(t),«(t)) + ||T(t)||i,M > W(V«(t)) - (J2?(t),tt(t)) + \\l(t)\\ 1)M , 

which implies ||7(£)||i,Ar > ||7(0lli,M- As ^(t) < 7(2) ^f ri_1 -a.e. on M, we deduce 
that 7(i) = 7(t) J^ n_1 -a.e. on M for every t G [0, T], which concludes the proof. □ 



Theorem 13 . 71 can be used to explain the mechanical meaning of the internal variable 
7 in the model case ip(x, y) := \y\. Indeed, if t > (u(t),j(t)) is a quasistatic evolution 
with 7(0) = and ip(x,y) := \y\, then ()3.4|) shows that j(t)(x) coincides with the 
maximum modulus of the amplitude of the opening reached by the crack at x up to 
time t. 

Remark 3.8. As 1 1— > 7(t) satisfies (J2.7|) by Theorem 13 .71 the mechanical interpreta- 
tion given in Section El shows that the term \\j(t) — j(0)\\i t M m f)3.3|) represents the 
energy dissipated in the crack process in the time interval [0, £]. 

Remark 3.9. In our model, the dissipation term in the energy functional comes from 
the expression ||7 V0([f]) — 7||i,m and is nonlinear in 7. This turns out to be the main 
mathematical difference between our model and the model considered by Mielke and 
Mainik and Mielke in ^21 Section 4.2] and Section 6.2], where the dissipation 
term is linear. 

We are now in a position to state our main result. 

Theorem 3.10. Let (u ,7q) G W 1,p (fi \ M;R m ) x L 1 (M) + be globally stable at time 
t = 0. Then there exists an irreversible quasistatic evolution t \— > (u(t),j(t)) such 
taa*(u(0) j7 (0)) = («o,7o). 



We introduce a notion of convergence for the functions 7, which is the counterpart 
of the notion of convergence of sets introduced in j^j. The main property of this 
convergence is that, if Uk converges weakly in W 1,P (Q \ M; M m ) to some function u 
and </>(M) < 7fc ^" _1 -a.e. on M, then <f>([u]) < 7 Jf n_1 -a.e. on M. 

Definition 4.1. Let 7^,7 G L°(M) + . We say that 7^ o^-converges to 7 if the 
following two conditions are satisfied: 

(a) if uj ->> u weakly in W 1 *^ \ M; R m ) and 0(N) < 7fej jT^-a.e. on M for 
some sequence kj — > 00, then 0([w]) < 7 ^ n_1 -a.e. on M; 

(b) there exist a sequence w l G W 1,p (Vt \ M; R m ), with sup^ </>([«*]) = 7 n_1 -a.e. 
on M, and, for every i, a sequence G W 1,p (£l \ M;M m ), converging to u l 
weakly in W 1 *^ \ M; M m ) as fc -> 00, such that 0(KD < 7fc ^ n_1 -a.e. on 
M for every i and /c. 

Notice that we do not require any upper bound in L 1 (M) + for the functions 7&. 



> 




^(t)(«(t),7(*))-^(0)(«(0),7(0))> 
(9W(V«(s)),V^(s)) - (Jzf(s),V>(s)) - (jS?(s),u(*)>) 



4. Some tools 



Remark 4.2. If 7^ a^-converges to 7, then in particular there are functions u\ and 
u l in W 1,p (Vt \ M; R m ) such that condition (b) in Definition 14. II holds. We define for 
every k 

7 4 := sup <KM) and 7^ := sup (f)([u{]). 

j=l,...,i J=l,— 1* 

With this notation it turns out that 

7 = lim 7* and 7^ > sup 7^,, 

for every k. 

Remark 4.3. If 7^ cr^-converges to 7, then 



7 < limsup7fc, Jtf? n -a.e. on M, 

as we can see by modifing the proof of Lemma 14.41 below. Notice that the inequality 
can be strict, even when 7^ converges pointwise to a function 7. As an example, 
consider n = 2, m = 1, p = 2, Q = ]-2, 2[ 2 and M = [0, 1] x {0}. Let lk G L°(M) + 
be defined as follows: 

7 fc (x) := < for i = 0, . . . , k - 1. 

[0 for *e 

It follows from homogenization theory (see [Zj, ^Hj, [IHj) that condition (a) in Defini- 
tional] is satisfied with 7 = 0, hence 7^ a^-converges to 0. Furthermore 7^ converge 
in measure to 1, so up to a subsequence we have pointwise convergence to 1 =: 7 > 7. 

We prove in the following lemma that the L 1 -norm is lower semicontinuous with 
respect to a^-convergence. 

Lemma 4.4. Let 7^, 7 G L°(M) + . If jk -converges to 7 then 

\\l\\i,M < liminf ||7fc||i,Af- (4.1) 

fc— >oo 

Proof. From the hypothesis it follows in particular that there are functions u\ and 
u l in W l,p (VL \ M; M m ) which satisfy condition (b) in Definition 14.11 With notation 
from Remark 14.21 let us prove that for every i 

||7i|i,M < liminf ||7^||i iM . (4.2) 



k— >oo 



Extracting a subsequence we may assume that liminf fc I^U^at is a limit. As [u 3 k ] — > 
[tt J ] strongly in L P (M; W 71 ) for j = l,...,z, we can extract a further subsequence 
such that [i^,] — > [iz J °] pointwise ^ n_1 -a.e. on M for j = 1, . . . , i. By the lower 
semi continuity assumption (<^ 3 ) this implies 

Y < liminf 7* Jf^-a.e. on M. 

k^oo 

By the Fatou lemma we obtain (|4.2jl . which yields 

IMIi.m < liminf ||7fc||i,Af- 

fe— >oo 

We then pass to the limit as i tends to infinity and obtain (|4.1jl . □ 

We now prove a compactness result for the notion of o~ v - convergence. 
Lemma 4.5. Every sequence in L°(M) + has a a v -convergent subsequence. 
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Proof. Let us denote the L p -norm by || • || p . Let 7 fc G L°(M)+, let w h G L°°(fi\M; R m ) 
be dense in L P (Q \ M; R m ), and, for every positive integers /, /i, and k, let us consider 
the problem 

mm{\\Vu\\ p p + e\\u-w h \\ p p }, (4.3) 

where the minimum is taken over all functions u G W 1,P (Q \ M; lR m ) such that 
<KM) < 7fe ^ n_1 -a.e. on M. 

To prove that the minimum is achieved, we take a minimizing sequence and we 
easily obtain that it is bounded in W ljp (Q \ M; M m ). Then, up to a subsequence, we 
can pass to the limit and by using our lower semicontinuity assumption (</? 3 ) we can 
prove that the limit function is actually a solution to the minimum problem (j4.Hj) . 
This solution, which is unique by strict convexity, will be denoted by . Notice 
that this function is bounded in W l,p (Q \ M; M m ) uniformly with respect to k, thus, 
up to a subsequence, we can pass to the limit in k and get that there is a function 
u t,h such that w fc ' — ^ -u £ ' h weakly in \ M; IR m ). Further we define 

7 := sup JT^-a.e. on M. (4.4) 

In this way point (b) of Definition 14.11 is automatically satisfied. 

We need to prove point (a). To this aim, let Vj — ^ v weakly in W l,p (Q \ M; M m ) 
be such that </>([fj]) < 7a^ ^ n_1 -a.e. on M for some sequence kj — > oo. We want to 
prove that <j>([v]) < 7 n_1 -a.e. on M. By density there is a subsequence of u>/j, say 
w^, which converges strongly to v in L P (Q \ M;M m ). Let £j — > +00 be such that 
— Whi\\p — >■ as i tends to infinity. By the minimality of u£ % we have 

l|V«Jj h< ie + £,||«Sf i -«; fci |g < IIV^^ + ^IK- -uU* 

Then w fc *.' ' is bounded in W /1,p (fi \ M; M m ) uniformly with respect to j, and passing 
to the limit as j tends to infinity we get 

|| W^H* + U\\u 1 ^ - w hi f v < sup ||V^-||£ + ii\\v - w hi \\ p p . 

Since £i\\v — w^Wp — > as i tends to infinity, this inequality ensures that Vu li,hl is 

bounded in L P (Q \ M; M. m ) uniformly with respect to i, and u£. ' — — > strongly 

in L p (fi s M;M m ). As w hi -> v strongly in L p (tt \ M;M m ), we deduce that w^ 1 
converges weakly to v in W 1,p {Vt \ M;M m ). Then converges strongly to [i>] 

in L P (M; W 71 ). Passing to a subsequence, we may also obtain pointwise convergence 
Jf? n_1 -a.e. on M. By flOJ) we have (j){[u k ' hi ]) < 7 JP^-a.e. on M, so that the lower 
semicontinuity assumption (^3) yields 0([f]) < 7 ^f 7ra ~ 1 -a.e. on M, which is precisely 
the conclusion to point (a) in the definition of a^-convergence. □ 

We shall use the following Helly-type compactness result. We recall that a function 
t 1 — > 7(t) from [0, T] into L°(M) + is said to be increasing if 7(5) < j(t) J# yn ~ 1 -a..e. on 
M, whenever < s < t < T. 

Lemma 4.6. Let t 1— > jk(t) be a sequence of increasing functions from [0, T] into 
L°(M) + . Then there exist a subsequence j^., independent oft, and an increasing 
function t 1— > j(t) from [0, T] into L°(M) + , such that 7^.(2) o~ p -converges to j(t) for 
every t G [0, T] . 

Proof. Let fbea countable dense subset of [0, T] containing and T. By Lemma I4~H| 
using a diagonal argument, we can extract a subsequence, still named 7fc(i), and an 
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increasing function t t— > y(t) from D into L°(M) + , such that 7fc(i) cr^-converges to 
7 (£) for every t E D. 
Let us define 

7(t+) := inf 7(5) and j(t— ) := sup 7(s), 

s>t,s£D s<t,seD 

for every t G [0, T\. It is easy to prove that: 

(1) 1 (t-) = ^(t) = 7 (t+) for every i G D; 

(2) 7(t-) < 7(t+) for every t G [0,T]; 

(3) if s < t, then 7(5+) < y(t-). 

Define £? := {t G [0,T] : 7 (t+) = 7 (i-) JT^-a.e. in M} and 7(f) := 7(t-) = 7 (t+) 
for every t £ E. Note that by (1) D is contained in E and the definition of j(t) 
agrees with the original one on D. Then the definition of a v - convergence and the 
monotonicity condition imply that 7 fc(i) a^-converges to j(t) for every t G E. 

Let us show now that the set E c := [0, T] \E is at most countable. For every pair of 
positive integers i, k we set A iyk '■— {t G [0, T] : \\ (j{t+)Ak) — (j(t—)Ak)\\i t M > ^/i}, so 
that we have E c is the union of the sets A^. Therefore it is enough to show that each 
set A itk is finite. Let t\ < • ■ ■ < t r G A ijk . Since, by (3), (y(tj-x+) A k) < ("f{tj—) A k) 
for j = 2, . . . , r, we get 

r 

T -<J2 U(7(*i+) A fc ) - (7(*i~) A < ||T(*r+) A fc||i,M < kJe n -\M), 

1 3=1 

so that r < ikJrff n ~ l (M), which implies that A^k is finite. It follows that E c is at 
most countable, thus we can conclude the proof of the lemma by applying again the 
compactness Lemma [4.51 for every t G E c , together with a diagonal argument. □ 

The following result plays a crucial role in the proof of point (a) in the Definition 13. 41 
of quasistatic evolution. 

Lemma 4.7. Let 7^, 7 G L°(M) + . Assume that 7^ cr^- converges to 7. Then for any 
v G W l ' p (Q \ M;R m ) with <f)([v\) G L 1 (M) + the following inequality holds true: 

hmsup \\{4>{{v\) -7*) + IIi,m < ||(0(H) -7) + I|i,m. (4.5) 

fc^oo 

Proof. It is not restrictive to assume that the limsup is a limit. Let u l and u\ be the 
functions considered in point (b) of Definition 14. 11 During the proof we shall use the 
notation introduced in Remark 14.21 As y\ < 7^ J$? n ~ 1 -a,.e. on M, we have 

(0(M)-7fc) + <(0(H)-7 fc ) + , 

hence 

lim ||(0(H) - 7fc ) + Hi,M < Bp inf ||(0(H) -7D + IIi,m. (4.6) 

fe^oo fc— >oo 

Passing to a subsequence, we may assume that [u\] converges to [u % ] M ?n ~ 1 -a,.e. on 
M. By the lower semicontinuity assumption (<^ 3 ) we obtain 

7* < lim inf 7* JT^-a.e. on M, 

fc^oo 

so that Fatou Lemma gives 

limsup ||(0(H) -7l) + ||i,M < ||(0(H) 
fc— +00 
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which, together with ()4.6|) . yields 

lim ||(0(M)- 7fc ) + l|i,M< ||(0(H)-Y) + IIi,m. 

fe— >oo 

As Y — > 7 J^ n_1 -a.e. on M, inequality ()4.5|) can be obtained by passing to the limit 
as i — > oo. □ 

Remark 4.8. The conclusion of Lemma [4. 71 does not hold, in general, when 7^,7 G 
L°°(M) + and 7^ — ^ 7 weakly* in L°°(M). Consider, for instance, the case n = 2, 
m = 1, Q — ]— 4,4[ 2 , M = [— 71, 7r] x {0}, and define 7 fc(V) := 1 + sin(fcxi), where X\ 
denotes the first coordinate of x. Then, 7^ converges to 7(2?) := 1 weakly* in L°°(M), 
but f)4.5|) is not satisfied for <f>([v\) = 1, since in this case ||(0([i/|) — 7fc) + ||i,M = 2 for 
every k, while - 7 ) + I|i,m = 0. 



5. The discrete-time problems and proof of the main result 

In this section we prove Theorem 13 . 1 ( )l by a discrete-time approximation. We fix a 
sequence of subdivisions (t\.)o<i<k of the interval [0,T], with 

= t°<4<---<^- 1 <^ = T, (5.1) 
lim max (4 - tj. -1 ) = 0. (5.2) 

k ~^oo l<i<k 

For i = 1, . . . , k we set Jg* = JSf(4), ^ = ^(f fc ), ^ = 

For every k G N we define and 7^. by induction as follows. Let («o>7o) be a 
minimum energy configuration at time t — 0. We set («°,7°) := («o>7o) an d define 
(u\, jI) as a solution of the minimum problem 

min{4!( M , 7 ):7eL 1 (M) + , 7 > 7 r 1 , u <E AD^)} , (5.3) 
where the inequality means that 7 > YfT 1 J? 11 - 1 -a.e. on M. 
Remark 5.1. Consider the minimum problem 

min {w(Vu) - (S*,u) + ||0(M) V ■ « = ^ on «9 ft}, (5.4) 

where m is assumed to belong to W 1,P (Q\M; M m ). Then the following two conditions 
are equivalent: 

(a) the pair (i4>7fc) is a solution to (j5.3jl : 

(b) 4 is a solution to (jED and 7^ := ^j. -1 V </>(KD ^ n-1 -a.e. on M. 

The existence of a solution of ()5.3|) (or equivalently 1)5.4)1 ) can be easily obtained by 
using the direct methods of the calculus of variations. The compactness of a mini- 
mizing sequence follows from ()2.3j) and positiveness of ip. The lower semicontinuity 
follows from (Wi), (W2), (^3), and from the compactness of the trace operator. 

For every t G [0, T] we define 

T k (t) = 4 u k (t) = u\, 7 fc (t) = 7|, Vfc(0 = ^(**)> 

where i is the greatest integer such that ^ < t. Note that Uk(t) = Uk(rk{t)), 7fc(t) = 
7k(n(t)), Mt) = ^t(t) = ^f(r fc (t)) and <f fc (t) = S{r k {t)). 
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Remark 5.2. Since ip k G ADfyl,^ 1 ), then by Remark we deduce that the 
LP- norms ||V^|| P and || u l k ||p are bounded uniformly with respect to i and fc. Passing 
to the piecewise constant functions t \— > Vufc(t) and t i— > Uk(t), we have that there 
exists a positive constant C such that 

||Vu*(f)|| p <C and \\u k {t)\\ p <C (5.6) 

for every fc and for every t G [0,T]. Since S'kit) (uk(t) , 1k(t)) is bounded uniformly 
with respect to k, we get also that 

||7fc(*)l|i,M<C, (5.7) 
for every k and for every t G [0, T]. 

We introduce now a sequence of functions which play an important role in our 
estimates. For a.e. t G [0, T] we set 

e k (t) := (0W(Vu fc (f)), V#)) - (JSf*(*),^(t)) - (J&?(*), «*(*)). (5.8) 

In the following lemma we present the main energy estimate for the discrete process. 

Lemma 5.3. There exists a sequence R k — > swcfc t/iat 

^fa(*))M*),7k(*)) < «r(0)(uo,7o) + / 9 k (s)ds + R k , 

Jo 

for every k and for every t G [0, T] . 

Proof. We need to prove that there exists a sequence R k — > such that 

4X4,71) <<nO)K 7 o)+ / 0fc(s)ds + i2 fc , 

Jo 

for any and for any % = 1, . . . , fc. 

Let us fix j and fc with 1 < j < fc. Since u-^T 1 = i\? h ~ on cfofi, and [i4 _1 + 
ipj, — ipjT ] = [u 3 k ~ ] ,J^f n ~ 1 -a.e. on M, the function u k ~ + ip k — ip 3 k ~ belongs to 

AD(-0£,7£~ ), hence S k { u k^k) - ^k( u l~ + ^1 ~ ^k~ >7iT )• The proof now can be 
concluded arguing as in the proof of [U Lemma 6.1]. □ 

We are now in a position to prove our main result. 

Proof of Theorem YS.lOX Let (t k ), < i < fc, be a sequence of subdivisions of the inter- 
val [0, T] satisfying (|5.1|) and (J5.2|) . For any fc consider the pairs {u k ,Y k ) inductively 
defined as solutions of the discrete problems (|5.3|) for i = 1, . . . , fc with the initial 
condition (m°,7°) = (u ,j ). Let Tk(t), Uk(t), 7fc(t), and ^(i) be defined by ()5.5|) for 
any t G [0, T]. By Lemma there exists a subsequence of 7fc(t), independent of t, 
which cr^-converges to ^ooif) G L°(M) + , for every t G [0,T]. By ()5.7j) and Lemma FQ1 
we have 7oo (t) G L X (M)+. 

Let 6k{t) be defined by ()5.8|) for a.e. t and let 

6>oo(t) := limsup0 fc (t). 

fc— >oo 

By (Q and (Q we deduce that 
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where || • ||* is the norm in the dual space of W l,p (Q \ M; W 11 ). Since the right-hand 
side of previous formula belongs to -^([O, T]), we deduce that 9^ belongs to i^ 1 ([0, T]), 
too, and using the Fatou lemma we get 

r T k(t) rt 
limsup / 6k(s)ds< I 9 00 {s)ds. (5.9) 



k- 

For a.e. t G [0, T] we can extract a subsequence 9^. of 9k, depending on t, such that 

floo (*) = lim M*)" 

By fl5.fij) the sequence ttfe.(t) is bounded in W 1,P (Q\M; MJ 71 ), therefore we can extract 
a further subsequence, still denoted by «&■(£), which converges weakly in W 1,P {VL \ 
Af; M m ) to a function ««,(*). 

Since </>([«&,; (£)]) < 7fe,-(*) ^ n_1 -a.e. on M, by point (a) in Definition 14.11 we 
have <fi([uoo(t)]) < l<x{t) =^ n-1 -a.e. on M. On the other hand, as itfc.(t) = ipkj(t) 
J^ n_1 -a.e. on <9 f2, we have also Uoo(t) = tp{t) J^ n_1 -a.e. on <9 ^, so that Uodt) G 

AD(V>(0,7oo(0) for ever y te[0,T}. 

The next step is to prove that the pair (Moo(^), 7oo(^)) satisfies property (a) of Defi- 
nition l3~4l To this aim, let 7 G L 1 (M) + , 7 > 7oo(£) and t> G AD(ip(t), 7). By the min- 
imality of the incremental solutions 7fc(0)> we have that (wk(£), 7fc(i)) < 
^t(*)(ffe, 7fc(*) V0([u])), where ffc := v+tjj^t) — if)(t). Since the functional w 1— > W(V-u) 
is weakly lower semicontinuous and strongly continuous, and the function t 1— > jSf(t) 
is continuous, we immediately obtain 

W(Vuoq(*)) < liminf W(Vu fc (t)), W(Vv) = lim W(Vv fc ), (5.10) 

fc^oo fc^oo 

(&(t), Uoo (t)) = lim (Jg? fc (t),u fc (*)), (J2f(t),«> = lim (JSf fc (t),t;*). (5.11) 

k^oo k— >oo 

So far we have easily obtained that 

W(Vttoo(t))- (jS?(t),«oo(t)) < 

< W(Vv) - (Sf(t),v) + limsup ||(<KM) -lk(t)) + \\i,Mi (5.12) 

fc^oo 

where the last term in right-hand side comes from the equality 

( 7 V0(H))-7=(0(M)-7) + , (5-13) 

which holds for every 7 G L°(M) + . In order to obtain that the pair (u^t), 7oo(£)) sat- 
isfies point (a) in Definition 13. 41 of quasistatic evolution we want to apply Lemma \4. 71 
To this aim we need to know that 0(["Uoo(^)]) £ L 1 (M) + . By ()5.7j) in Remark 15.21 
we have that ||7fc(t)||i,M is bounded uniformly with respect to k. As Ukit) belong 
to AD(ipk{t), lk(t)), the sequence <fi([uk{t)]) is bounded in L 1 (M) + , and by the lower 
semi continuity assumption (<^ 3 ) we obtain that (j>([uoo(t)]) G L 1 (M) + thanks to the 
Fatou lemma. Then we can apply Lemma 14.71 and we get 

W(VUoc(*))- (&(t), Uoo (t)) < 

< W(Vv) - (Jf(t),v) + ||(0(H) -7oo(t)) + ||i,M. (5.14) 

Applying (|5.13|) to the last term in the right-hand side of (|5.14|) we conclude that 
^(*)(«oo(*),7oo(f)) < £(t)(v, lo c(t) V0(H)) < 7) for every t G [0,7] and 

point (a) of Definition 13.41 is satisfied. 

By the definition of the discrete problems, for every k the function t 1— > 7fc(t) is 
increasing. Passing to the a^-limit, the same property holds for t 1— > 7oo(0> so that 
point (b) of Definition 13.41 is satisfied. 
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It remains to prove point (c). For a.e. t define 

9(t) := (dW(V Uoo (t)),Vij>(t)) - (J?(t),m) - (Sf{t), Uoo (t)). 

Arguing as in the proof of Theorem 3.15] we get 

M*) = 0(t), (5.15) 

for a.e. t G [0, T\. This in particular means that the map t \— > 6(t) is measurable. 
Since we have proved that for every t G [0,T] the pair (Uoo(^), 7oo(^)) satisfies points 
(a) and (b) of Definition 13 A\ we are in a position to apply Theorem 13.61 and get 

^)Ko(t),7co(t))-^(0)K),7o) > / 0{s)ds. 

Jo 

By (@~H), f|5~TU|l . and dSHJ) we get 

<r(t)(uoo(*),7oo(*)) < liminf ^(t)K(t),7fc/t)) < limsup S k {t){u k {t), 7fc (*)). (5.16) 

J-* 00 fc^oo 

Using Lemma [5.31 and taking ()5.9|) and (J5.15j) into account, we obtain 

limsup 4(t)(« fe (t), 7fc (t)) < <f(0)(u ,7o) + / 0(s) ds. (5.17) 

fc^oo Jo 

By flSjH and (j5~T7jl we get that 

^(*)(uoo(*),7oo(*))<^(0)(txo,7o)+ / 0{s)ds 

Jo 

holds true for any t G [0,T], and this concludes the proof. □ 

In the following theorem we prove that for every t G [0, T] the energy for the 
discrete-time problems converges to the energy for the continuous-time problem. We 
emphasize that the theorem is true for any irreversible quasistatic evolution t i— > 
(tt(t),7(£)) corresponding to a given t i— > 7 (t), not only for the one obtained as limit 
of the solutions of the discrete-time problems. 

Theorem 5.4. For every t G [0, T] /et «&(£) and 7jfe(i) be defined as in the beginning 
of the proof of Theorem \H.l(K Assume that jk{t) o%- converges to 7 (t) G L 1 (M) + for 
any t G [0,T]. Let t i— ► (-u(t), 7 (t)) be an irreversible quasistatic evolution. For a.e. 
t G [0, T] Zei 9k(t) be defined as in and set 

e{t) : = (dw(vu(t)),vm) - (j?(t),m) - (jgf(t), «(*)>. 

T/ien 

W(Vu(t)) - (jSf(t),ti(t)> = lim (W(Vu k (t)) - (JZ> k (t),u k (t))), (5.18) 

fc— >oo 

|| 7 (*)l|l,M = lim ||7fc(*)l|l,Mj 
fc— +00 

for every t G [0, T\. Furthermore 

e k ^e m^dcT]), 

so that there exists a subsequence of 9 k which converges to 6 a.e. in [0,T]. 
Proof. For the proof we need to show that 

lim W(Vu k .(t)) = W(V Uoo (t)), (5.19) 

for every t G [0, T], where u kj (t) is the subsequence constructed in the proof of 
Theorem 13.101 and u^t) is its limit. To this aim, let Vj := u^it) + ip kj (t) — ip{t). 
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By the minimality of the pair (itfc,,.(0> 7%(0) we obtain that <^-(0( M fcj(0> 7fcj(0) — 
<^k (t)( v ji lkj(t) V (j)([uoo(t)])) , and passing to the limit as j goes to infinity, we get by 
(P , and (I5TTD 

limsup [W{Vu kj {t)) - (Sf k .(t),u kj (t))] < 



< limsup [W(V^) - (J^.(t),%> + 11(0(^(0]) -7^(0) + IIi,m] = (5.20) 

SUp ||(0([Woo(O]) — 7fc j (0) + l|l,M- 

Since 7fc J -(t) u^-converges to 7oo(t), ^y Lemma E~TI we have 

limsup ||(0([uoo(t)]) -7fc j (0) + lli,M < 0. (5.21) 

Taking into account (|5.2()|) and (|5.21jl we get in particular that 

limsupV^Vu^f)) < W(Vmoc(0). (5.22) 
jr'^oo 

This, together with (f5TT0"j) . gives (jSHJ). 

To conclude the proof it is sufficient to follow the arguments of the proof of [3J 
Theorem 8.1]. □ 

The result can be improved under strict convexity assumption. 



Theorem 5.5. In addition to the hypotheses of Theorem \5.4\ assume that £ 1— > 
W(x,£) is strictly convex for a.e. x G Q \ M and i/ioi y l— ► <p( x ,y) ^ s convex for 
,ff n - x -a.e. x G M. Then u k {t) -> it(t) strongly in W^ P {Q \ M;lR m ) ; /or every 
t G [0,T]. 

Proof. We observe that for every t G [0,T] and 7 G L X (M) + the functional v 1— > 
<^(t)(i>,7) is strictly convex on the set of functions t> G W /1,p (f2 \ M; M m ) with t> = 
■0(0 ^ n_1 -a.e. on <9ofi. Therefore for every t there exists a unique function u G 
AD (ip(t), 7(0) suc h that the pair (it, 7(i)) is globally stable at time £. It follows that 
u(t) coincides with the function Woo(0 constructed in the proof of Theorem l.ll ( )l and 
that the whole sequence u k (t) converge to u(t) weakly in W 1,p {Vl\M; W 71 ). Therefore 
flTTHjl implies that W(V« fe (t)) -> W(Vu(t)). Using [TEl Theorem 3] we deduce that 
Vu k (t) — > Vit(t) in measure. As 

I Vu fc (i) - Vw(t)| p < 2 p - 1 a 1 [W(Vu k (t)) + W(Vu(t))} + 2 p - 1 a 1 b , 

the conclusion follows from the generalized dominated convergence theorem. □ 

6. EULER CONDITIONS 

In this section we study the Euler conditions satisfied by globally stable pairs 
(it, 7) G W lj, (Q s M;M m ) x L X (M) + . Let us fix t G [0,T] and let (it, 7) G W^fl \ 
M; R m ) x L X (M)+ be globally stable at time t, and let v G W lj> (Q \ M; M m ) be such 
that d = ^ n-1 -a.e. on 9qO. Hence for every e > the function it + eu belongs to 
ADijpit)^ V 0([it] + e[u])), and by the global stability of the pair (u, 7) at time t, we 
have that S{t){v,,i) < £{t){u + ev^ V </>([it] + e[v])), therefore 

liminf ^(t)(n + ei;,7V0(N+eM))-^(t)(n,7) > g (g 1} 

e^0+ e 

The weak formulation of the Euler conditions will be obtained from this inequality. 
Without loss of generality, we assume that jSf(t) is given by (|2.5|) . and we omit the 
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dependence on time. After some standard calculation, one can express ()6.1j) in the 
following form 



/ (deW(x, V«) - H) : Vv dx - [ fvdx- I gvd^- 1 ^ 
- [ W + 9 e v e ) + liminf m[u] + £[V]) ~ 7)+lkA/ > 0, 



(6.2) 



for any v G W x *(p, \ M; M m ) such that v = Jf^-a.e. on <9 n. 

To continue our analysis we need now to specify the form of the function ip. More 
precisely, we consider ip: M x M. m — > [0, +oo] defined by 

<p(x, y) := <po(x) + <p{%i y) for y ^ and <p(a;, 0) := for all a; G M, (6.3) 

where <£>o € L 1 (M) + and p: M x W 71 — > [0, +oo] is a Borel function. We assume that 
for every x G M the following properties hold: 

(1) <p(x, y) — if and only if y = 0; 

(2) the function £(ar, ■) belongs to the space C°(M m ) n C^R" 1 \ {0}); 

(3) p(x,0) = 0; 

(4) there exists an L°°-function such that \d y <p(x,y)\ < <p(x) for any y ^ 0, 
where d y p(x, y) denotes the vector of the partial derivatives of <p with respect 
to y; 

(5) the limit 

ip(x,y) := \im d y ip(x,ey)y (6.4) 
exists and is finite for any y ^ 0. 
Remark 6.1. By using de l'Hopital Theorem, one obtain immediately that 

ib(x,y) = lim ^ ' , 

for any x G M, ?/ ^ 0. It follows from the positiveness of <p that ip > 0. Moreover, 
we get easily that ip is positively 1-homogeneous with respect to y, i.e., ip(x,Xy) = 
\ip(x,y), for every A > 0. Furthermore, by ()6.4|) and (4), we get also 

\i>{%,y)\ — ^(^OM f° r every x G M and y 7^ 0. (6.5) 

The main result of this section is a theorem which makes explicit the Euler con- 
ditions obtained from ()6.2j) in the case of the function ip specified above. Before 
stating the theorem, we establish a general result concerning closed linear subspaces 
of L*(f2), for an arbitrary Radon measure /x on Q. We will apply this result to the 
measure \i = ,ffl n ~ x L M. 

The characteristic function of any set E is denoted by 1^, i.e., Ie{%) = I ii x E E, 
1e(x) = otherwise. 

Lemma 6.2. Let /j be a Radon measure in Q and let Y be a closed linear subspace 
of L^(Q) with the following properties: 

(a) if u, v G Y , then 11V v &Y ; 

(b) ifueY and oj G C™(Q), then uou G Y. 

Then there exists a Borel set E c fi such that Y = {u G L^Q) : u = fi-a.e. on E}. 
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Proof. We begin by proving that 

if u G L (ft) and \u\ < \v\ for some v G Y, then u G Y. (6.6) 

Indeed in this case there exists u G L£°(f2) such that u = uv and there is a sequence 
cjfe G C^°(Q) such that ^ is bounded in L™(Q) and — *■ u /i-a.e. on Q. By (b) we 
have u k v G F, and by the Lebesgue dominated convergence theorem ui k v — > ojv = u 
in L Xl Xi). Since F is closed, we conclude that u G Y. 
Now we prove that 

if u G y and i > 0, then u A t G y and (u - t) + G K (6.7) 

As |w A t\ < \u\, we have u A i G y by (j6.6|) . Since (w — t) + = u — u A t, we obtain 
that (w - £) + G y. 
Next we prove that 

if u G y and t > 0, then l{ n>t } G y (6.8) 

where {u > t} := {x E Q : u(x) > t}. By fl6.7j) we deduce that for every > we 
have k(u— t) + Al G Y. As [£;(«-£)+] Al — > l{ u >t} pointwise and [k(u—t) + ] Al < |w|/t, 
the convergence takes place in _L^(f2) and we conclude that l{ u> t} g y. 

Let (uk) be a sequence dense in Y and let E be the intersection of the sets {u k = 0}. 
It is easy to prove by approximation that u = /i-a.e. on for every u £ Y. 
Conversely, let u G L^(Q) with w = /i-a.e. on £7. For every let 

4:={miVm 2 V---Vmjc> 

By (a) and ()6.8|1 we have 1^ G y, so that (M^ ) A u + and (ZcIaJ A u~ belong to Y, 
by (jfi.fij) . As (felAft) A m + — >■ m + and A w~ — >■ w~ in L^(O) we conclude that 

m g y. ' □ 

Lemma 6.3. Let D C M, let Yff" be the set of all functions of the form [v], with 
v G W 1 ' p (tt\M; R m ) and [v] = Jf^-a.e. on D, and letY% be the closure ofY£ in 
L 1 (M n dM; R m ). Then there exists a Borel set D (unique up to M' n ~ 1 -equivalence), 
containing D, such that Y% = {w G L\M \ dM; R m ) : w = Jp-'-ce. on £>}. 

Proof. Let be the set corresponding to the case m = 1. It is easy to see that 
Y^ 1 = (Yd) 771 . Therefore it suffices to prove the lemma in the case m — 1. 

The conclusion follows from Lemma 16.21 applied to Yd . It is enough to verify that 
conditions (a) and (b) are satisfied by Yd. Condition (b) is trivial. To prove (a) we 
consider an open set U Cfl \ M, with C l boundary and M C dU, such that U lies 
on the negative side of M. Given two functions u and v G W 1,P (Q \ M) it is easy to 
check that [u] V [v ] = [u V (v — v + u)], where u and v G W 1,P (Q) coincide with u and 
v on {7, respectively. □ 

In the following theorem we will consider a function u G W 1,P (VL \ M; M m ) such 
that the divergence of the matrix field d^W(x, V«) — if belongs to L q (Q s M; M m ). 
It turns out that its normal trace (9^W(a;, V«) — if)z/ is defined as an element of 

(W l ~p' p (diVL] W 71 )) 1 '. Moreover, we have that the normal traces (d^W(x, Vu) — H)®v 
and (d^W (x, Vtt) — H) e u (defined on the positive and negative side of M) are 

both elements of the space (W p' p (M \ dM; R m ))'. The duality pairing between 

(W l ~v> p (M s dM; W 1 ))' and W 1 ~p' p (M \ dM; R m ) will be denoted by (•, •). 
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Theorem 6.4. Let t G [0, T] and (n, 7 ) G W 1,P {VL \ Af; M m ) x L 1 (M)+ 6e #/o&a% 
stable at time t. Assume that (p: Af x M m — > [0, +oo] zs defined as above in \6.S\) and 
it satisfies (l)-(5). Then 

-div(d^W(x, Vn) - ff ) = f on Q \ Af, (6.9) 
(d s W(x, Vm) - H)v = g on <9 x fi, (6.10) 
{dzW(x, Vn) - if)®!/ + g® = (9 £ W(z, Vm) - ff ) e z/ - o e on M \ dAf. (6.11) 



Let ns define 



A 
B 
D 



{x G Af : < </>(M)(z) = 7(x)}, 

{x G M : = 0([u])(x) and 7(2;) = <fo(x)}, 

{xeM : 7 (x) < <^o(^)}, 



and Zet .D 6e the set associated with D by Lemma 16791 Then there exists h G 
L°°(Af \ D; R m ) snc/i tfmi 

((%W(x,Vn) -H)®v + g®,[v]) = I h[v}dJ^ n '\ (6.12) 

for every v G W 1,p {Vt \ M; IR m ) snc/i £/ia£ [n] = J^ n_;L -a.e. on f7 Moreover 

(a) for Jrf? n ~ 1 -a.e. x G A\Z) i/ie vector h(x) belongs to the segment joining and 

(b) /or ^ n_1 -a.e. x £ B \ D the vector h(x) belongs to the bounded convex set 
K(x) := {a G R m : ay < ^(x,y), Vy G M m }; 

(c) /or je 71-1 -a. e. x G Af \ (A U 5 U D) we fcaue = 0. 

Remark 6.5. It is easy to see that, if D is (^""-"--equivalent to) a closed set, then 
D = D (up to ^"^-equivalence). A more difficult proof shows that the same result 
is true if D is (^""-"-equivalent to) a quasi closed set with respect to (l,p)-capacity. 
It is clear that, if ipo = 0, then D = D = 0. 



Remark 6.6. For Jf? n G Af the vector h(x), obtained in Theorem 

represents the cohesive force exerted from the positive lip of the crack on the negative 
lip. The theorem shows the conditions satisfied by the cohesive force on the different 
regions of Af determined by the respective relations between </>([it]), 7 and (po- 



Proof of Theorem 6^7 Since 0([n]) < 7 J4? n 1 -a.e. on Af, we have (</>([«]) — 7) + = 



J^^-a.e. on Af. If [v] = Jf n_1 -a.e. on Af, then the liminf in (Q is actually a 
limit and it is zero. Therefore ([6.9)1 . (|6.10jl . and (|6.11|) can be obtained from (|6.2|) 
by standard argument involving integration by parts and a suitable choice of the test 
function v G W 1,p (Q]M. m ). 

To shorten the notation, we set h := (d^W(x, Vn) — H)®v + g® on Af \ 9Af . As 

explained before the statement of the theorem, we have h G (W ~p' p (M \dM ); K m ))'. 
So far, we may rewrite (|6.2j) as 



(-ft, M) + liminf nv^M+ e H)-7) + lli,M > Q> (613) 

e-»0+ £ 

for any v G W 1)P (fi \ Af; R m ) such that v = J^^-a.e. on c? fi. 
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Let us extend the definition of ip by setting ip(x, 0) = for every x G M. Now we 
prove that 

[u] +€W) -7) + l|l,M 



lim 

£ (6.14) 
{d y <p{x, [u])w) + l A + 4>{x, w) 1 B ) dJ^ n -\ 

I M 

for every w G L l {M \ dM;R m ) with w = Jf n-1 -a.e. on D. To this aim, it is 
convenient to split the set M into the union of the following two disjoint subsets 
A' := {x G M : [u](x) ^ 0} and B' := {x G M : = 0}. 

On A', as 0([m]) < 7 ^ n_1 -a.e. on M, we have that 

(0(M+ew)-7)+ < (0(M -0(M))+ = (y(x, M + ew) - <p{x, [u}))+ < 

£ ~~ £ £ 

< (<^(a;)w) + , 

J^ n_1 -a.e. on M, where we used (jfi.3j) . and assumptions (3) and (4). Moreover, we 
have that 

(0(N+£ g W) " 7)+ - (d y 0(x, [u])w)+ 1 A J^-a.e. on A', 

because A = {0 < = 7}. By the Lebesgue dominated convergence theorem we 

get 

as £ -> 0+, for every w G L l (M \ <9M; M m ). 

Let us consider now the integral over B' . If u> G L X [M \ <9M; M m ) and to = 
J$f n ~ 1 -&.e. on D, we have 

(<p(ew) - 7)+ _j 
= U J?r -a.e. on JJ, 

e 

thus we can focus on the set B' \ D. As 7 > tpo Ji? n ~ 1 -&.e. on M \ D, for every 
w G L X {M \ <9M; R m ) with w = Jf^-a.e. on D, we obtain 

(0(ew) - 7) + < {(j>{ew) - ip ) + = <p(x,ew) < 

J^ n_1 -a.e. on M \ D, where we used (lfi.3|) . and assumptions (3) and (4). Moreover, 
by Remark|H3 we get that 

^ £W )~^ + 1b = ^ w )i B jf"-i. a . e . on B', 

as £ -> + , for every iu G L X (M \ <9M;IR m ) with w = Jf^-a.e. on D. We can 
apply again the Lebesgue dominated convergence theorem and obtain 



J B' £ JM 



as e -> 0+, for every it; G L l {M \ <9M; lR m ) with iu = Jf n_1 -a.e. on D. This con- 
cludes the proof of ()6.14|) . We note that this equality cannot be true if the condition 
w = J4? n ~ 1 -a,.e. on D is violated, because in this case 

||0(£w) -7) + l|i,Af r HOo + <p(ew) -7) + II x,M 
lim = hm = +00. 

e^o+ e e^o+ e 
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Let be the space defined in bemmalOl Notice that C W l ~p' p (M\.dM); R m ). 
By (jnHj) and (l6~Tl) we have 

(-h,w)+ I [(d y <p{x, [u])w) + l A + tP(x,w)l B ]djr n - 1 > 0, (6.16) 

J M 

for any w G Yq\ In order to localize this inequality, we prove first ()6.12|) . Due to our 
assumption (4) and to ()6.5|) . if we apply ()6.16|) to w and — w we deduce that 

\(h,w)\ < \\(f\\oo\\w\\i,M^D , (6.17) 

for every w G Y™. It follows that there exists a function h G L°°(M \ .DjIR" 1 ) such 
that 

for every u> G Yp\ This implies that ()6.12|) is satisfied. By density from ()6.16|) we 
obtain 

dje n ~ l > 0, (6.18) 

for every w G Kg 1 . Since by Lemma I6~31 we have Kg 1 = {w G L l (M \ <9M; IR m ) : 
w; = ^ n_1 -a.e. on /}}, we conclude that 

-h(x)y + (d y 0{x, [u](x))y) + l A {x) + #(x, y) l B (x) > 0, (6.19) 

for every y G W 1 and for n_1 -a.e. x £ M \ D. 

In particular, for ^ n_1 -a.e. x £ A\D the equality d y (p(x, [u](x))y = implies that 
h(x)y = (it is enough to use (|6.19|) with y and — y), so that for a given x G A \ I) 
the two vectors d y ip(x, [u](x)) and /i(x) are parallel, hence there exists X(x) such that 

h(x) = \(x) d y p(x, [u](x)) for JT^-a.e. x £A\D, (6.20) 

and it is easy to verify that < A(x) < 1, by using again ()6.19|) . In this way we get 
condition (a). 

On B \ D, from ()6.19j) we obtain 

-h(x)y + #(x, y) > for JT^-a.e. x G B \ D, (6.21) 

for every y G K m , which is precisely condition (b), by the definition of K. On the 
remaining part of M \ D, from ()6.19jl we get condition (c). This concludes the 
proof. □ 

Remark 6.7. If (fo(x) > for Jf n_1 -a.e. x G M, and (u, j) = (u(t),j(t)) for an 
irreversible quasistatic evolution, then ()3.4|) implies that the set B \ D is nonempty 
only if there exists y G M m \ {0} such that <p(x, y) = 0, for some x G M. This 
happens, for instance, in the Griffith model, where <p is given by ()2.6p with a > and 
6 = 0. In this special case, condition (b) becomes h(x) = ^f n ~ l -a,.e. on B \ D, 
because K(x) = {0}. 

Remark 6.8. If for every x the functions £ t— > iy(x,£) and y ip(x,y) are convex, 
then for any t G [0, T] and 7 G L 1 (M) + , the functional u 1— ► ^(t)(u,7 V </?([«])) is 
convex. Therefore, it is possible to prove by standard arguments that conditions (a), 
(b), and (c) of Theorem 16.41 are equivalent to the inequality 

- I hwd^ + lim HWM+^)-7) + lli,M 
Jm £ ^ 0+ e 



hw + (d y ip(x, [u])w) + 1a + tp(x, w) 1 



B 
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for every w G Y%. Thus, Euler conditions ()6.9|) . ()6.10p . (j6.11|) . (a), (b), (c) are not 
only necessary, but also sufficient to global stability. 

We show now an example of a scalar problem, where the Euler conditions of The- 
orem 16.41 lead to a simplified set of boundary conditions. 



Example 6.9. Let m = 1, p = 2, W(x,g) := ||£| 2 , H{t) := 0, g®(t) = g G {t) 



0. 



4>(y) '■= \y\, which correspond to the energy functional: 
1 
2 



<f(i)(u, 7 ) 



\Vufdx+ / 'ydJtf' 

Q\M J M 



re— 1 



f(t)udx- [ g{t)udJ^ n - 1 . 
J din 



Let t G [0, T] and (u, 7) G H^' 2 (Q \ M) x L 1 (M) + be globally stable at time t. Then 
we are in a position to apply Theorem 16 .41 and the final part of Remark 16.51 obtaining 

(-Au = f(t) onfixM, 
u = %fj(t) on doQ, 

on diQ, 

on MH {0 < |[«]| < 7}, 
< 1 and ^[w] > on M n {|[u]| = 7}. 



9m . . 
du 



dv 



du , 



9z/ 



By Remark 16.81 we have also that if u solves the previous boundary value problem for 
a given 7, then the pair (u, 7) is globally stable at time t. 



7. The case of linear elasticity 

In this section we show that, with some modifications, it is possible to consider also 
the case where the uncracked part of the body is linearly elastic, which is excluded 
by the first inequality in (|2.1jl . 

Let p = 2 and m — n > 1. We assume now that the bulk energy relative to the 
displacement u G W 1,2 (Q \ M; M. n ) has the form of linear elasticity 

A{x)Eu : Eu dx, 

where Eu := ;(Vu+ (Vw) T ) is the symmetric part of the gradient of u, and A satisfies 
the following properties: 

(Ei) for every i6(], A(x) is a linear symmetric operator from the space 

symmetric n x n matrices into itself, and the map x 1— > A(ar) is measurable 
(£2) there are two positive constants Cq and C\ such that 



nxn _r 
sym ^' 



coiei 2 <^)e^<ciiei 2 



(7.1) 



irnxn 
sym " 



for every a; G \ M and £ G 

For the sake of simplicity in the notation we introduce the C 1 map Q: L 2 (Q \ 
M; M^) -> R defined by 

Q(*) := / : 
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for every f 6L 2 (fi\M; M££), whose differential dQ : L 2 {tt \ M; M%%) -> L 2 (fi \ 
M;M^) is given by 

(<9Q(#),$) = 2 / dz, 

for every $, $6L 2 ((]\ M; M™^), where (•, •) denotes now the scalar product in the 
space L 2 (fix M;M££). 

For every t G [0, T] the total energy of an admissible configuration (it, 7) G iy 1,2 (fi\ 
M,M n ) x L 1 (M)+ at time t is now defined as 

<f(t)(«,7) := Q(£u) " + hlli.M- 

Once we have the energy functional, we introduce the notion of global stability as in 
Definition 13.11 

Since the (n — l)-dimension of <9ofi is positive, Korn inequality holds (see, e.g., 
|Hj): there exists a constant C = C(fi, <9ofi) such that 

||Vu|| 2 < C\\Eu\\ 2 for all u G W l ' 2 {Vt] W) such that w = 0on <9 fi. 

As an immediate consequence, we get the following Korn- type inequality: 

||Vtt|| a < C\\Eu\\ 2 + (C + 1)|| W>|| 2 (7.2) 

for every u G W 1 ' 2 ^ \ M; M n ), and ip G iy 1,2 (fi; M n ) such that u = ip on <9 fi. 

Thanks to (|7.2|) . we still have an a priori bound for the displacement u as in 
Remark 13.21 

The definition of irreversible quasistatic evolution of minimum energy configura- 
tions is now given replacing (dW(Vu(t)), Vip(t)) by (dQ(Eu(t)),Eip(t)) in Defini- 
tional 

Thanks to the Korn-type inequality ()7.2j) . Theorems 13. 7\ 13.101 and 15.51 (and 
Remark l3.5|) continue to hold, with essentially the same proofs, if we replace W(Vu(t)) 
and (dW(Vu(t)),Vip(t)) by Q(Eu(t)) and (dQ(Eu(t)),Eij)(t)), respectively, and a 
similar substitution is done for Uk{t). 
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